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We have determined the electron-coupUng spectrum of superconducting Bi2Sr2CaCu2 08+(5 from 
high-resolution angle-resolved photoemission spectra by two deconvolution-free robust methods. As 
hole concentration decreases, the coupling spectral weight at low energies <15 meV shows a twofold 
and nearly band-independent enhancement, while that around ^-^65 meV increases moderately, and 
that in ^130 meV decreases leading to a crossover of dominant coupling excitation between them. 
Our results suggest the competition among multiple screening effects, and provide important clues 
to the source of sufficiently strong low-energy coupling. Ale ~ 1, in an underdoped system. 
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The coupling of electrons with other excitations plays 
an essential part in possible pairing mechanisms of su- 
perconductivity, and it concomitantly makes an electron 
appear as a slower and heavier quasiparticle. In reality, 
the electron is coupled with multiple kinds of excitations 
of various frequencies. Hence, the energy resolved data 
on the electron coupling provide important clues to the 
pairing glue. For high-Tc cuprates, it is believed that the 
strong electron correlation comes from the proximity to 
the Mott insulating phase i2j. However, the behav- 
ior of the group velocity of renormalized quasiparticle is 
intriguing. With decreasing hole concentration, the de- 
crease in velocity does not occur on the energy scale of 
the Mott transition [SHS], but at far lower energies, <40 
meV, as reported recently [6l [7] . The mechanism of this 
nontrivial mass-enhancement remains far from clear. In 
order to pin down the source of the coupling strength and 
the pairing attraction, we have to unravel the multilevel 
renormalization effects. Therefore, a thorough investi- 
gation of the energy, doping and band dependences of 
electron coupling is required. 

Angle-resolved photoemission spectroscopy (ARPES) 
is an excellent tool for studying the interaction from 
the electron side [7HT4]. Extracting the coupling infor- 
mation from subtle features in ARPES data has been 
attempted by the maximum- entropy method [8^ |9j and 
Richardson-Lucy method ^3 However, such de- 

convolution is a severe integral inversion problem, and 
possibly sensitive to statistical noise [15 . Here, we have 
developed two deconvolution-free robust methods for de- 
termining the "effective" coupling spectrum from high- 
resolution low-temperature ARPES spectra, noting the 
causal nature of the mass-enhancement factor defined as 
A(co') = —{d/duj)J^{uj)^ where T^{uj) is electron self-energy 

m- 



In this paper, we report a systematic low-energy 
ARPES study of the electron-coupling spectrum of super- 
conducting bilayer cuprates, Bi2Sr2CaCu208+5. Quanti- 
fying the impacts of three coupling features on the quasi- 
particle mass, we show that their contrasting evolutions 
with hole concentration cause a change in the dominant 
coupling excitation to occur. We propose possible sce- 
narios for the mass enhancement with underdoping. 

High-quality single crystals of Bi2Sr2CaCu208+5 were 
prepared by the traveling-solvent floating-zone method 
and a post annealing procedure. The hole concentration 
p has been deduced from Tc using a phenomenological re- 
lation, Tc/T^^'' = l-82.6(p-0.16)^ where T^^"" = 91 K 
[17 . Hereafter we label the samples by the doping level, 
i.e., underdoped (UD), optimally-doped (OP) or over- 
doped (OD), combined with the value of Tc. The ARPES 
spectra were collected at BL-9A of the Hiroshima Syn- 
chrotron Radiation Center using a Scienta R4000 electron 
analyzer. Instrumental energy and momentum resolution 
was 5 meV and 0.004 A~^. The samples were cleaved in 
situ^ and kept under an ultrahigh vacuum (pressure under 
5 X 10~^^ Torr) at T = 9 K during the measurements. 

Figures [TJa)-(d) show the low-energy region of ARPES 
spectra. Despite a difficulty in controlling surface qual- 
ity, a tiny nodal bilayer splitting was resolved even for 
the UD samples, as demonstrated in Figs. [TJe) and (f). 
The full momentum width at half maximum, 0.009 A-i, 
for UD66 is narrower than the previous studies [71 113] . 
We found that the spectral intensity ratio between the 
bilayer-split bands drastically changes with photon en- 
ergy hv. We adopted hv = 8.1 and 7.0 eV for simul- 
taneous observation of the bonding band (BB) and the 
antibonding band (AB) in Figs.[TJa)-(c), and for selective 
observation of AB in Fig. [TJd), respectively. 

The quasiparticle group velocity is given by the slope 
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FIG. 1: (a)-(d) Energy- momentum plots of ARPES spectra 
along the nodal direction of Bi2Sr2CaCu208+(5 for UD66 (un- 
derdoped, Tc = 66 K), 0P91 (optimally-doped, Tc = 91 K), 
and OD80 (overdoped, Tc = 80 K) samples, (e) Momentum 
distribution curves (MDCs) at cj = for hv — 8.1 eV. (f) 
Energy distribution curves (EDCs) at /c = k^^ for hv ^ %.\ 
eV. (g) Quasiparticle dispersions determined by MDC fitting 
for UD66 (red), 0P91 (purple), and OD80 (blue). Open and 
filled circles denote the results for hv — 7.0 and 8.1 eV, re- 
spectively, (f) Inverse Fermi velocities of BB and AB, l/vp^ 
(diamonds) and l/t'p^ (triangles), and momentum-to-energy 
ratio of bilayer splitting width, A:bs/<^bs (filled squares). 



of dispersion, v^{{jj) = duj/dk. As hole concentration de- 
creases, the splitting narrows more rapidly in energy than 
in momentum as shown in Figs.[T]^e),[TJf), andjljh), pro- 
viding a clear evidence for the decrease in Fermi velocity. 
The narrowing tendency with underdoping is consistent 
with what is expected from the decrease in out-of-plane 
conductivity [18]. In Figs. [TJa) -(d), a dispersion kink is 
consistently observed at \uo\ ~ 65 meV [5,^5^. Moreover, 
the dispersions at low energies, \uj\ < 15 meV, are sub- 
stantially curved for UD66, whereas they seem relatively 
straight for 0P91 and OD80. We can rule out the effect 
of transition-matrix elements, because the dispersions de- 
termined with different photon energies are identical as 
shown in Fig. [ijg). Comparing the Fermi velocities of 
BB and AB in Fig. [ijh), we find that BB exhibits an 
effective-mass enhancement similar to AB with under- 
doping [7|. These results indicate that the low-energy 



interaction is nearly independent of the bilayer bands. 

Thus, we have determined the quasiparticle dispersion 
and scattering rate over a wide energy range by imposing 
cj-independent bilayer-splitting parameters on the fitting 
analysis of momentum distribution curves (MDCs). Fig- 
ures [2| a) and|2|e) show the peak position k{uj) and natu- 
ral half- width A A: (a;), respectively. Approximating bare- 
electron velocity vq by a constant, we obtain the forms, 
k{uj) = [uj — ReT>{uj)]/vo and Ak{uj) = —lmE{uj)/vo. Fig- 
ure Isja) shows that a small dispersion kink at ~10 meV 
[inillll El evolves with underdoping into a large and 
broad feature around ~15 meV. The difference between 
the inverse group velocities at the Fermi level and at 40 
meV shown in Fig. [2|d) indicates that the coupling with 
low-energy excitations is abruptly enhanced upon enter- 
ing the underdoped region. 

The real and imaginary parts of mass enhancement, 
l-\-X{uj) = Vo/vg{uj), are directly deduced from the energy 
derivatives of k{uj) and A/c(a;), respectively, by the forms, 

dk{uj) _ 1 + ReA(cj) _ 1 



duj 
dAk{uj) 
duo 



ImA(cj) 

^0 



Vg(cj)' 



The differential scattering rate, ImA(a;), at T = repre- 
sents a kind of coupling spectral function, which includes 
the effects of k and u dependences of electronic spectral 
function. Note that the uj dependence of 1 + A (a;) is ir- 
relevant to the uncertainty of '^o unlike that of Rei;(co'). 

The energy dependence of mass enhancement ReA(u;) 
is presented as l/vg{uj) in Fig. |2j^b). A steplike mass 
increase at ~ 65 meV is the typical behavior of electron 
coupling with a certain bosonic mode, and in good agree- 
ment with the results of optical studies [19]. By contrast, 
a cusplike mass enhancement at a; = for UD66 appears 
to have no saturation of slope, as shown in Fig. |3|b). 
This indicates that the onset of coupling spectral weight 
is quite close to = 0, and may cause virtually singular 
behavior of the quasiparticles. 

The effective coupling spectra ImA(co') have been de- 
duced from the Kramers-Kronig transform (KKT) of 
{d/duj)k{uj) and directly from {d/duj)Ak{uj)^ as shown 
in Figs. [2|c) and |2|f), respectively [20^. Furthermore, 
dividing the integral. 



A = ReA(O) 



f 

Jo 



ImA(a;) 



duj, 



into three energy parts, we have quantified the partial 
coupling constants as shown in Figs. [2|d) and|2|g). Al- 
though the experimental accuracy is better for the peak 
position k{uj) than for the peak width A A: (a;), the deriva- 
tion is more direct in Fig.j2|f), and the slope of the spec- 
tral background in Fig. I^Tc) remains uncertain due to 
the extrapolation required for KKT ^]. To this ex- 
tent, the results from k{uj) and Ak{uj) are consistent. As 



3 




FIG. 2: (a) Quasiparticle dispersion determined from MDC-peak position, k(uj) = [cj — Kel2{uj)]/vo. (b) Inverse group velocity, 
l/vg(uj) = [1 + KeX{uj)]/vo, determined from [d / duj)k{uj) . (c) Effective coupling spectra, —1vsi\{ijj) / vq ^ deduced from the 
Kramers-Kronig transform (KKT) of —{d/duj)k{uj) [20 . (d) Doping dependence of inverse group velocities at Fermi level, 1/vf 
(black circles), and at 40 meV, 1/vao (black triangles), and partial coupling constants, X/vq (color), deduced from the KKT of 
{d / duj)k{uj) . (e) Quasiparticle scattering rate determined from MDC-peak natural half width excluding instrumental resolution, 
/S.k{uj) — —lmT^{ijj)/vQ. (f) Effective coupling spectrum, — ImA(aj)/vo, directly determined from —(d/duj)Ak{uj). (g) Doping 
dependence of partial coupling constants, X/vq, deduced from {d/duj)Ak{(jj). (h) Strong electron-phonon coupling at low hole 
concentration, (i) Weak electron-phonon coupling at high hole concentration. The low-energy (LE, red circles), intermediate- 
energy (IE, blue triangles), and high-energy (HE, green diamonds) parts are defined as 4 < < 40 meV, 40 < \uo\ < 130 
meV, and 130 < |cj| < 250 meV, respectively. The differential coefficient at uj has been evaluated within an energy window of 
cj — W{uj) < cj < cj + W{uj), where W{uj) = 5.5 + 0.15|c(;| meV, from simple difference between both ends of the window (color 
dots) and by least-squares linear regression method (black curves). The data are offset for clarity in (a), (b) and (e). 



hole concentration decreases, whereas the intermediate- 
energy (IE) part around the 65-meV peak shows moder- 
ate increase jSmUS], the low-energy (LE) part abruptly 
shows twofold enhancement with a curve quite similar 
to that of l/v-p. This behavior is manifest not only in 
the curvature of the dispersion k{uj) in Figs. [2|a) and 
[sja), but also in the slope of the scattering rate Ak{uj) in 
Fig.[2|e), and is consistent with the cj-linear term of scat- 
tering rate deduced from tunneling spectra [21 . Further- 
more, Figs.[2|d) and|2|g) show that the impact of the LE 
part exceeds that of the IE part on the underdoped side, 
implying a crossover of the dominant coupling excitations 
from ~65 meV to <15 meV. Assuming ^'o ^ 4 eVA based 
on the local-density-approximation calculation [22 , suffi- 
ciently strong coupling of order unity, A^^ ~ 1, is realized 
only with the low-energy excitations. 

The energy-dependent enhancement of the coupling 
spectrum with underdoping is likely related to high- 
energy electron-electron interaction [9 . In Fig. [2|[f), one 
finds that the substantial coupling weight linear in uj 



extends over \uj\ > 130 meV beyond the phonon cut- 
off. Moreover, Fig. [2|g) suggests that the high-energy 
(HE) part decreases to zero towards the superconductor- 
to-insulator transition point. Therefore, this part should 
be ascribed to the electron-electron interaction expressed 
as ImS®~® oc cj^. With sufficient hole concentration, 
the bare Coulomb potential is quickly screened by these 
high-frequency electronic excitations, and thereby the 
retarded response of low-frequency excitations is sup- 
pressed. Such suppression becomes more drastic as the 
frequency of the coupling excitation decreases. Thus, 
the contrasting behaviors of the LE, IE, and HE parts 
are consistently interpreted as the competing effect in 
screening the Coulomb potential. 

A candidate for the origin of the large effective mass in 
the underdoped system is the bare strong coupling with 
low- frequency optical phonons as illustrated in Fig. |2|h) . 
Rameau et al. assigned the small low-energy kink of 
optimally-doped Bi2Sr2CaCu2 08+5 to the c-axis phonon 
involving the vibration of Bi and Sr atoms [I0l[23] . At the 
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FIG. 3: Dispersion anomaly at low energies, (a) Energy devi- 
ation from the straight line which intersects with the exper- 
imental dispersion at a; = —40 and meV. Filled and open 
marks denote the result for hv = 8.1 and 7.0 eV, respectively, 
(b) Inverse group velocity l/vg{(jj) for UD66, compared with 
a logarithmic function, 0.151n|c<;| (dotted curve). 



node of d-wave gap, the low-energy scattering is only pos- 
sible without in-plane momentum transfer. The phonons 
of out-of-plane momentum can provide such forward scat- 
terings, and thus are compatible with the (i-wave super- 
conductivity also as a pairing interaction [24 . 

In addition, the effective mass may also be enhanced 
near the quantum critical point between the metallic 
and insulating phases, as proposed from the quantum 
oscillations [6]. The instability toward some competing 
orders of charge or spin induces the divergent critical 
fluctuations [25H27] , even though the optical phonons 
make no contribution to the mass divergence. Phe- 
nomenologically, the LE parts of ReA(co') and ImA(co') for 
UD66 resemble the marginal- Fermi-liquid form, X{uj) cx 
In IcJc/^l — 1 — ^7r/2 for < where oOc is cutoff energy, 
as shown in Figs. (sfb), |2];b), |2];c) and||;f) [25 . 

Finally, we note the effect of elastic forward impurity 
scattering. Tunneling spectroscopy measurements have 
revealed that nanoscale areas, where the coherence peak 
is absent, emerges for p < 0.13 [28l |29]. The local de- 
pletion of quasiparticles implies breakdown of the static 
screening, thus giving rise to the static-potential inho- 
mogeneity that originates from out-of-plane distant im- 
purities [29 , and would result in a dramatic increase in 
elastic forward scattering. Given the conelike dispersion 
around the node of the d-wave gap, the elastic scattering 
rate may have quasilinear energy dependence [HI [131 EO] • 
In that case, the energy and mass of nodal quasiparticle 
may also be renormalized by the second-order forward 
scattering process. 

In conclusion, our ARPES study has revealed the 
nearly band-independent and contrastingly energy- 
dependent evolution of the electron-coupling spec- 
trum with hole concentration. In underdoped 
Bi2Sr2CaCu208+5, the strong coupling weight (Ale ^ 1) 
is distributed around ~15 meV with an onset at a; 
0. As hole concentration decreases, the LE part shows 



twofold rapid enhancement, the IE part increases mod- 
erately, and, by contrast, the HE part decreases to zero 
towards the superconductor-to-insulator transition point. 
This behavior suggests the competition among multiple 
screening effects as a possible origin of the mass enhance- 
ment. In terms of the dominant coupling excitation, a 
crossover from ~65 meV to <15 meV occurs upon un- 
derdoping, whereas the electron-electron part in ^130 
meV increase in presence with overdoping. The balance 
among these multiple interactions provides a new per- 
spective on the phase diagram of cuprates. In particu- 
lar, whether the low-energy interaction is pair-breaking 
or pair-binding would be an important subject of future 
study. 
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